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Abstract 



We study how the DGP (Dvah-Gabadadze-Porrati) brane affects particle 
dynamics in linearized approximation. We find that once the particle is re- 
moved from the brane it is repelled to the bulk. Assuming that the cutoff for 
gravitational interaction is M^, ~ 1/e, we calculate the classical self energy 
of a particle as the function of its position. Since the particle wants to go 
to the region where its self energy is lower, it is repelled from the brane to 
the bulk where it gains its 5D self energy. Cases when mass of the particle 
m < Svr^M^, and m > Stt'^M^ are qualitatively different, and in later case one 
has to take into account effects of strong gravity. In both cases the particle 
is repelled from the brane. For m < Stt'^M^ we obtain the same result from 
the 'electrostatic' analog of the theory. In that language mass (charge) in the 
bulk induces charge distribution on the brane which shields the other side of 
the brane and provides repulsive force. The DGP brane acts as a conducting 
plane in electrostatics (keeping in mind that in gravity different charges repel). 
The repulsive nature of the brane requires a certain localization mechanism. 
When the particle overcomes the localizing potential it rapidly moves to the 
bulk. Particles of mass m > Svr^M* form a black hole within 1/M* distance 
from the brane. 



-mail: mk679@nyu.edu 



1 



I. INTRODUCTION 



A phenomenologically acceptable five dimensional brane world theory with one infinite 
extra dimension was recently developed in The low scale of quantum gravity M^, is 

pulled (renormalized) at the brane by the high scale Msm that describes the brane localized 
standard model. More precisely, the 4D Einstein Hilbert term, with strength ~ M| ~ M'smi 
is induced on the brane. This effect ensures that the observer on the brane sees weak 4D 
gravity (Newton constant Gat ~ 1/Mp) up to the distance = Mp/M^. At distances bigger 
than Tc gravity becomes five dimensional. At short distances gravity is modified by quantum 
corrections at M~^. Short distance gravity measurements exclude modification of laws of 
gravity at distances bigger than ~0.1mm~ 1/lO^^eV Cosmological observations on the 
other hand suggest that gravity is not changed to distances of order ~ lO^^mm. Thus, 
present knowledge about gravity constrains the scale of gravity in this class of models to the 
range 

IQ-^V < M, < lOMeV. (1.1) 

Relativistic corrections, and the question how are they encoded in the tensor structure 
of the graviton propagator were studied in . The cosmological consequences of the model 
and especially the fact that the model gives rise to the accelerated universe (as observed 0) 
were considered in 0. The relevance for the solution to the cosmological constant problem 
was considered in 

In the present paper we study how the induced 4D Einstein Hilbert term affects dynamics 
of the particle of mass m at distance yo from the brane. This question is important in order 
to identify experimental signatures of collider black hole production in this class of models. 
Since the particle itself induces the metric, there is no static metric in which one could study 
the behavior of geodesic lines (that describe the motion of test particle). Instead of trying to 
solve the problem in full relativistic theory, we will limit ourself to Newtonian approximation 
and use some basic facts about black holes. The result that we find is that the brane repels 
particles into the bulk where they have lower (bigger in magnitude and negative) self energy. 

In the next section we derive and briefiy discuss Newtonian potential. In the third 
section we find the dependence of self energy of a particle as a function of its distance from 
the brane. Once particle leaves the brane, the gradient of self energy forces it to go from 
the brane to the bulk. We also describe the process of black hole formation for particles 
with mass m > Svr^M^,. In the fourth section we present the 'electrostatic' analog derivation 
of repulsive force for particles with mass m < Sir'^M^. Finally, in discussion, we address 
questions regarding phenomenological consequences of the repulsive nature of the brane. 

II. NEWTONIAN POTENTIAL 

Action for the model is the sum of the 5D Einstein Hilbert term and induced the 4D 
term on the brane 

S = Ml j d^xdyVcn^^^ + M^i f d^x^/\g\R. (2.1) 
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Here we divide 5D coordinates into 4D part (Greek indices) and extra coordinate y like 
= [x^^ y), Gab is 5D metric and 7^(5) its curvature scalar, g^vi.x'^) = Gab{x^, U = 0)5^5^ 
is induced 4D metric (we take straight brane) and R corresponding scalar curvature. The 
tension of the brane is taken to be zero. If we take the limit of slowly varying weak fields, 
equations of motion reduce to the equations for deviation of the (700 component from flat 
space constant value (scalar gravity). The equation for (Euclidean) Green's function for the 
scalar gravity case reads 

(□4(1 + r,5{y)) - dl) G{x - xo, y, yo) = S\x - xo)S{y - yo), (2.2) 
which has the solution (see [Q]) 

Gip,y,yo) = -e~P\y-y'^\ - -e-^^M+M) ] . (2.3) 

p p 1 + 1 / ''^cP 

Let us evaluate the exact Newtonian potential at the point (x, y) due to a static source of 
mass m located at the position (x', yo). The potential is given as a Fourier transform of the 
Green's function ( p.3| ) integrated over the time 

in / I 1 

V{r, y, yo) = 



2rr, 



167r2M3 + {y - y^)^ + {\y\ + \yo\)^ 

' (l^l+l^ol-)A'= (Tom + \yo\-^r)/r,)~e''^'/'^rom + \yo\ + ^r)/r,))\ . (2.4) 



Here r = \x — x'\ and ro(-2) is an incomplete gamma function (see appendix). The potential 
can be expanded in powers of l/vc 

1 m 1 r 

^ = oTT? — arctan— j — r (2.5) 

Wir^M^rrc \y\ + ml 



V- = i M±M arctan + 1 ('A in + ^1^1 + 



2 

IQir'^MlrTc \ \y\ + \yo\ 2 Vrc/ 

(2.6) 



where 7 ~ 0.5772 is Euler's constant and superscripts on potential denote terms in expan- 
sion. The potential to first order in l/r^ was discussed in detail in 0. Let us briefly discuss 
potential ( [2.41 ). If the mass is on the brane [yo = 0), the first two terms in ( ^.4[ ) cancel. 
The potential on the brane {r,y = 0), at distances r <C r^, is four dimensional and the 
Newton's constant is G = l/(327rM^rc). As r/vc increases towards one, the second term 
in (|2.6|) weaken its strength. Finally, for r ^ Vc {pvc ^ 1) the potential becomes purely five 
dimensional (the first term in equations ( |2.3| ),( P~^ )). Similar behavior occurs if one looks 
at the potential at (r = 0,y). For small ?/ it is a weak four dimensional potential with 
constant 2G/ti. For y/rc^l one can find the form of the potential by expanding ( p.4|) 
(see appendix) and again obtain the expected five dimensional behavior. Up to a constant, 
potentials have the following short distance expansion and asymptotic behavior 



3 



in 1 / 2 7^ T \ Tfl 

V(2/ = 0,r«r,) = -— -— — 1 + In- , V{y = ^ r,) = (2.7) 

V{r = 0, 1/ « r,) = - 1 + ^ In A , = q, y » r,) = 

(2.8) 

Similar expansions can be easily obtained for a potential at any 'angle' in the r — y plane. 

If the mass is in the bulk we have two different cases. For particles on opposite sides 
of the brane, the first two terms in (|2.4| ) cancel and particles interact via the weak four 

dimensional gravity at distances ^(|y| + |?/o|)^ + ^ ^c- That means that the brane is 
shielding one side of the brane from the five dimensional gravitation of sources on the other 
side of the brane. The effective radius of shielding is ~ rc. If the sources are on the same side 
of the brane, the interaction is dominated by the first two terms in ( p.4|) . Masses, sufficiently 
far from brane, basically interact via strong five dimensional gravity. 

One can illustrate this behavior by plotting the contours of the constant potential of 
the body as it moves from the bulk towards the brane (Fig.l). At Schwarzschild radius 
(700 ~ (1 + 21^) diverges. Although we don't have a relativistic solution to the system, one 
would expect that the Schwarzschild surfaces (black hole horizons) behave as surfaces of 
constant potential V ~ —1/2. 




0.5 1 1.5 2 2.5 



Figure 1: Surfaces of constant potential V = — (l/167r^)m/M*. Point mass is on r = 0, yo = 
1, 0.5, 0.2, 10"^, 10"'^, 10~^ in units of M~^. Inset: anisotropic cutoff distance e for a mass on the brane 
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III. SELF ENERGY 



In this section we will evaluate the classical self energy of the particle of mass m in the 
presence of the DGP brane. Our main assumption, along the lines of [Q, is that the gravity 
is cut off at distances e ~ 1/M*. Classically, gravitational self energy is determined by the 
cutoff distance and the form of the potential. Since the potential changes with the position, 
the gravitational self energy of a particle will be a function of its distance from the brane. 
Gradient of self energy will give rise to a force that will try to move particle to the region 
of lowest gravitational self energy. For m < Svr^M* we will use Newtonian approximation, 
since the potential is weak at the cutoff distance ~ M^^. For m > Svr^M*, we will use a 
Newton-like approximation, which will incorporate some basic facts of general relativity. 



A. Case m < Svr^M* 



For m < Stt^M* Schwarzschild radius in the bulk r^s = {l/M^)^m/8Tc'^M^: ^ is smaller 
than the cutoff so Newtonian approximation is justified. Let us remind how we calculate 
gravitational self energy of a particle. One wants to find an energy gained by assembling a 
particle of mass m. Equivalently, one can find the energy needed to destroy a particle by 
taking away infitesimal pieces of matter and removing them to infinity, where the potential 
is defined to be zero. If the force starts to act at a distance e from the center of mass 
distribution of a particle, one then finds the expression for self energy to be 

W = - f\l-M)dM r^^^dr = -V{e). (3.1) 



Je dr 2 

In our normalization of Newton constant, that leads to self energies in four dimensional 
theory (on the brane) and five dimensional theory (infinitely far from the brane) 

If we take that the gravity cutoff e is just an inverse scale of gravity M^, then the ratio of 
self energies in pure 4D and 5D theories is 



2 VMp 



(3.3) 



In pure four or five dimensional theories energy of the mass density p{x) is given hy W = 
(1/2) / p{x)V{x)d^x. One can then use Gauss law to relate mass density and divergence of 
the field. After partial integration energy can be written as an integral over the square of 
the field ly ~ — / | W(x)p(i"x. By integrating energy stored in the field in pure 4D and 5D 
theories one again finds expressions (|3l^). We must stress here that the discussed theory is 
neither purely 4D or 5D theory. In particular, for the mass on the brane. Gauss law is not 
valid (if we call the 5D radial distance from the mass at uq = 0, then the field drops like 
~ l/r|, while the surface area increases like ~ r|). For this reason we will use ( ^.Ij ) when 
calculating self energy of the mass at an arbitrary position in 5D space (one should not use 
Gauss' law). 
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Let us look at the self energy of the particle on the brane. By using prescription ( p.l|) 
we find that the self energy is that of a four dimensional particle. However, since our 
space is not isotropic there is an apparent ambiguity in self energy. It depends on the 
direction from which we assembled particle. Let us define polar coordinates p = + and 
= arctan(r/?/) and assemble particle at ?/o = 0, r = by bringing infitesimal masses from 
infinity and direction from the bulk . Self energy of the particle on the brane then varies 
by a factor 7r/2 (same as Newton constant) for angles < < 7r/2: 

W(y, = 0) = -^J^f^' (3.4) 

SZtt^ \Mp/ sm0 e 

Since the self energy must be a well defined quantity, we conclude that the factor 0/ sin0 
defines the physical cutoff distance when the mass is on the brane. Certainly, the space 
in question is not isotropic, and we cannot assume that the gravity cutoff surface is a 3- 
sphere, but rather a surface of constant field strength with average distance from the particle 
~ (inset to Fig.l). For the arbitrary position of particle j/o, the expression for self 

energy to first order in l/vc is 

T../ X 1 fm'^ \ ( ^ (. 1 \ 1 1 ( esin0 

W{y,) = TTIT 1 - ^-T7-7-^T-T7-7-r-— + —--Cretan 



4(yo/e)^ + 4(?/o/e) cos0y r^esin^ \2?/o + ecos< 

(3.5) 

For yo > e, dependence on can be neglected and the self energy has the form 

W{yo) = wJl- } ] + W,^^ (3.6) 
V l + 4(We)V 7r(yo/e) 

For < yo < € (neglecting the small l/vc contribution of 4D self energy) , one gets the same 
answer by following argument. The first term in ( p.5|) represents isotropic 5D interaction 
and should be cut off at the surface of 3 sphere of radius e. The second term is anisotropic 
contribution and its cutoff has to be defined so that it doesn't depend on the angle 0. If 
we define (angle) dependent cutoff for the second term as = e'^f{yo/e,(f)) one finds that 
fiyo/e, 0) = (— 2(?/o/e) cos0+ \Jl + 4(?/o/e)^ cos^ 0). Plugging this back to (|3?5|) one obtains 
the behavior of 5D contribution as in (137 



To summarize, at the brane, the particle has 4D self energy, upon leaving the brane, 
within a couple of e distance, it gains the biggest part of its 5D self energy and looses its 4D 
self energy. The particle at yo > ^ will feel strong force 

dW{yo) 1 I ( yo\ 



dyo 647r2M3y3 

This force will try to push the particle to the bulk where its self energy increases in magnitude 
by the large factor of (Mp/M,)^. 

B. Case m > Svr^M* 

If the mass of the particle is bigger than the scale of gravity M^, we cannot calculate 
the self energy by cutting off the Newtonian potential at e ~ M~^. The reason is that 
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Schwarzschild radius in the bulk is bigger than the inverse scale M^,, and at distances shorter 
than Schwarzschild radius gravity is not weak. In following considerations we will not 
write negligible corrections of the order of l/r^. 

Let us calculate the self energy of a 5D black hole by approximating the black hole as 
an object that gravitates via Newton's law at distances bigger than event horizon r > rs5- 
What happens with potential at distances below event horizon doesn't influence the energy 
of the world outside horizon. Again we construct the self energy by assembling black hole 
at the origin out of infitezimal pieces of matter located at infinity. To assemble the point 
mass M^, we need to bring matter to a cutoff distance M~^, since the Schwarzschild radius 
is smaller then the inverse cutoff. This contribution equals to W^lm = M^) = — M^,/(327r^). 
Work gained in bringing the rest of the mass (from to m) would be the work done in 
bringing it only to the Schwarzschild radius. Outside the horizon there is no change in energy 
no matter how the potential changes below the horizon. Now we introduce a loose definition 
of the Schwarzschild radius as a radius at which Newtonian potential has the value —1/2. 
Since the potential energy gained by bringing the mass dM from infinity to Schwarzschild 
radius is by (our) definition dM/2, the total self energy of a 5D black hole of mass m is 

W{m > M„ — ^ oo) = W^im = M,) - - / dM = -^--^M, - - (m - M,) . (3.8) 

In our simplified model of a black hole, self energy is negative and of the order of (factor of 
1/2) the mass of a black hole. It is interesting to note that self energy could, in principle, 
be equal to the rest mass so that it wouldn't cost anything to produce it. The situation is 
reminiscent of the fact that the total mass of universe, Newton's constant and the Hubble 
radius conspire in such a way that it might not cost anything to create particles at the center 
of the universe, since their rest mass energy is of the order of their gravitational (negative) 
energy. 

Let us see how the formation of the black hole happens as we remove a particle of mass 
m > Svr^M^, from the brane. On the brane, the Schwarzschild radius rs4 ~ Mpi^(m/Mpi) is 
much smaller than the cutoff distance M~^. For that reason the self energy on the brane is 
given by W4 (|3.2| ). The self energy on the brane is much smaller than the self energy far away 
from the brane ( |3.8| ) and can be neglected. Thus, particles of mass m > Stt^M* will (as well 
as particles with m < Svr^M^,) be repelled from the brane to the bulk where their self energy 
is lower. However, the character of the repelling force will differ from the case of m < Svr^M^,. 
Let us define the Schwarzschild surface as a surface on which V{r,y,y(y) = —1/2. When the 
particle is removed from the brane, the Schwarzschild surface expands from a point (actually, 
the three sphere of radius rs4, ~ Mpi^(m/Mpi) ) anisotropically (Fig.l). After the particle 
reaches a certain value of yo, the Schwarzschild surface crosses the three sphere of radius 
that describes the gravity cutoff radius (Fig.2) . This crossing first happens for the 
value of = 0. Up to that point the self energy and the force on particle are the same as 
in the case m < Svr^M,,. After that point we cannot consider self energy cutoff distance 
and the self energy evolves different from the case m < Svr^M*. Moving the particle further 
into bulk, the Schwarzschild surface grows and takes over the M^, sphere at larger angles and 
finally, for some critical value of yo, the M^, sphere becomes completely contained inside the 
Schwarzschild surface. At that point we can say that the black hole formation is finished 
and the particle will have (up to corrections of the order of M^,) the self energy of a 5D black 
hole. 
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0.5 1 1.5 2 2.5 

Figure 2: The figure shows the formation of the black hole horizon in the r — y plane. The horizontal 
axis measures distance from the brane y and the vertical axis measures the distance along the brane r, both 
m units oiM-^. The mass of the particle is taken to be m = IGtt^M, i.e. twice the critical mass. Solid lines 
represent Schwarzschild surfaces, dashed lines represent surfaces of three spheres of unit radius and masses 
are represented by dots. The horizon emerges when the particle is at yQ~° = 0.207 and completely encloses 
the M, sphere one unit farther. 

Points where the Schwarzschild surface crosses the M* sphere for = 0, vr can be ex- 
pressed from potential ( p.4|) 

For m ^ Svr^M^, these expressions become ?/o~° ~ ^ix'^ jm and y^~^ ~ 1/M*. Thus, the 
horizon starts forming at ~ l-n"^ jm and is formed precisely farther. The force felt by 

the particle on < 2tx'^ /m is the same as in the case m < Svr^M^, {yoM^, ^ 1) 

Fy{yo < 21? jm) = -^m^M.yo. (3.10) 
For 2iT^/m < yo < 1/M*, neglecting terms of order in self energy, force is approximately 

F,(yo > 27rVm) ^ ^. (3.11) 

To summarize, for m > Sn'^M^, the 5D Schwarzschild radius is bigger than the cutoff 
M~^, and one cannot use the Newtonian theory to calculate the self energy of the particle. 
Modeling a black hole as an object that gravitates with Newtonian potential outside the 
horizon, we calculated self energy and estimated the force felt by the particle. As in the 
case m < Stt^M* particles are repelled from the brane. 
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IV. CONDUCTOR ANALOGY 



In this section we will rederive results of the previous section for particles with m < 
Stt^M^. from a different point of view. The lagrangian of our theory can be thought of as 
the lagrangian for a purely 5D theory with specific type of source localized aX y = 0. One 
would expect this correspondence to be valid everywhere except at y = (worldvolume of 
the source) because the source itself is a kinetic term for the 4D theory. At y = 0, the value 
of delta function diverges and the 4D kinetic term becomes dominant. Our gravity theory 
in this approach {y ^ 0) becomes equivalent to a 5D gravity in the presence of an infinite 
3-plane with a specific mass (charge) distribution. In the Newtonian approximation theory 
is equivalent to the electrostatic setup of a charge near the conducting plane. We will use 
symbol E for the gravitational field and sometimes interchange the terms mass and charge. 

Let us take mass m at position (r = 0,yo) and look at the field at position (r, y). One 
can ask what kind of charge distribution on the plane would produce potential ( p.4|) . 



r 



m 



y 



Figure 3: Gaussian pillbox for determination of the effective charge distribution p{r, yo) induced by 
the charge m at position (r = 0, yo) 



A component of the field in the y direction is discontinuous at y = with discontinuity 
(to a first order in I/tc) 



m \yo\- {r^ + yo)/{2r^) 



(r2 + 2/2) 



2\2 



(4.1) 



Applying the Gauss theorem on the 5D pillbox, as shown in Fig. 3, we can find the charge 
distribution on the plane 



VE 



4M3 



pix,y) 



Pir, yo) 



+ 



m 



1 m\yo\ 

7r2(r2 + yg)2 ' 27r2rc(r2 + yg) ' 



(4.2) 



The first term in the expression for charge density ( [4.2|) represents the distribution of negative 
charge, sharply peaked around r = 0, i.e. the projection of the charge position on the brane. 
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Its integral over the volume of the brane is independent of i/q and equals precisely — m. 
As Ho — > this term approaches distribution of a point-like charge — m5(r). The second 
term in ( [4.2| ) is a small (notice l/r^ suppression) distribution of positive charge, much less 
localized then the negative charge distribution. In order to find the total induced charge 
density we should integrate the exact expression obtained from the potential ( |2.4| ) that is 
correct to all orders in 1/rc. The integral of the charge distribution due to the third term 
in ( |2.4| ) can be obtained numerically and is equal to m. Thus, the total charge induced on 
the brane due to the charge m in the bulk is zero, as one would expect. 

Having the charge distribution, we can calculate the interaction energy between the 
mass and the induced distribution ('image' distribution of mass —m and the background 
distribution of mass +m). For y ^ our theory is just the 5D Newtonian gravity so the 
potential energy of interaction is (to first order in l/vc) 



W{yo) = 2 j p{x)V{x)dV = -^^^ J = IW (.MfJ ^ I' - T: 

(4.3) 

Cutoff effects in this derivation were neglected, so it is understood that yQ > e. The result 



( ^731) coincides with ( |3.6|) and gives the same force (|3.7|) . 

Let us summarize what happens in our 5D analog picture. Charge m in the bulk induces 
negative 'image' charge distribution of total charge — m, localized at the r = 0, and the small 
uniform background positive mass distribution. The total induced charge on the brane is 
zero. As charge m approaches the brane {yo — > 0) image charge distributions tends to the 
distribution of a point like charge — m, which strongly repels mass m. Finally, charge m 
and the image —m annihilate and m distributes itself uniformly on the brane. Described 
process is completely analogous to the behavior of the charge near the conducting plane. 
The only difference is that in electrostatics charges of the opposite sign attract and in gravity 
they repel. Thus, mass m in the bulk is repelled from the brane by its image —m. In this 
sense the DGP brane acts as a gravity conductor, shielding the fields and giving rise to a 
repulsive force. One could imagine constructing tensionless objects with this property that 
would gravitationally repel masses or act as gravitational dipoles. In cosmological setups 
the tensionless DGP brane would gravitationally shield ('shadow') parts of the universe and 
could modify cosmological evolution. 



V. DISCUSSION 

In previous sections we showed that the particles are repelled from the brane that induces 
kinetic term. By analogy with an ordinary wall with a tension [|T0| , we can loosely say that 
the induced kinetic term creates localized energy momentum density on the brain in which 
repulsive tension dominates over attractive energy density. Theories with low scale of gravity 
predict collider production of black holes. Because of the repulsive nature of the brane, black 
holes produced in collider experiments would be repelled to the bulk. 

The repulsive nature of the brane requires certain localization mechanism for standard 
model particles. We can distinguish two different cases. In the first case standard model 
particles are entities that cannot exist independent of the brane. Well known examples 
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are goldstone modes of broken translational invariance (elastic waves of the brane), modes 
of open strings with endpoints stuck on the brane, or simply fermionic zero modes on 
the soliton-like wall. In this case, particles feel force but they cannot escape to the bulk. 
Other possibility is that the standard model particles are entities that exist independent 
of the brane. Then, we have to introduce localizing potential AW that keeps them on 
the brane. Since on colliders we don't see events in which particles just disappear, the 
depth of the localizing potential AW would have to be bigger than ~ ITev. From the 
present bound on the size of universal extra dimensions one knows that the range of the 
localizing potential should be less than 300Gev |TT]. Localizing potential can be due to 



short range (contact) interactions with the matter of the first type, or the brane itself. For 
phenomenologically acceptable energy densities on the brane, the gravitational attraction 
can not provide localizing potential. The particle localized on the brane will feel effective 
potential which is a combination of short distance localizing potential and repulsive potential. 
With the potential of depth AW, all particles lighter than AW will be in stable equilibrium 
on the brane. Particles heavier than AW would be in a metastable state on the brane, 
because their self energy in the bulk is roughly their mass ( |3.8| ). Metastable particles can 
then tunnel through the barrier into the bulk. The brane can, in principle, be populated 
with both types of particles. Intrinsically brane particles would be stable (with respect to 
escape to the bulk decay). Particles trapped on the brane, on the other hand, can decay 
by escape to the bulk. Upon leaving the brane those particles would gain energy of the 
order of their mass (bulk self energy) in the vicinity of the brane (distance ~ 1/M*). The 
recoil effect of the brane would produce stable particles (goldstone modes, zero modes) with 
a total energy of the order of the mass of particle that escaped to the bulk. This kind of 
decay to the bulk would make missing energy signal on colliders smaller than one expected 
in a scenario with an ordinary brane. 

So far, all the discussion has referred to a brane of infitezimal thickness (delta function 
type brane). Real, physical branes, have finite thickness. It would be interesting to see how 
the finite thickness affects particle dynamics and if it can provide localization mechanism. 
Arguments that we used in the derivation of the repulsive force in 5D model, apply equally 
well to branes in space with more than one extra dimension. To completely understand 
particle dynamics, one would certainly like to have an exact relativistic solution. 
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APPENDIX A: INCOMPLETE GAMMA FUNCTION 

Incomplete gamma function Ta{z) is defined as 

fOO 

^a{z) = / e-H^-^dt. 

J z 

It satisfies 

^Tq{z) = - — , / To{z)dz = -e-' + zVo{z). 
az z J 

For small and large values of argument it has following expansions 

ro(^) ^ --i - \nz + z + 0{z^), z — ^ 
ro(^) ^ e-' (l/z - l/z^ + 2/z^ + 0{l/z^)) , z — > oo. 
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